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Abstract 

We investigate the localization of a hydrophobic - polar (HP) - regular copolymer at a selective 
solvent-solvent interface with emphasis on the impact of block length M on the copolymer behavior. 
The considerations are based on simple scaling arguments and use the mapping of the problem 
onto a homopolymer adsorption problem. The resulting scaling relations treat the gyration radius 
of the copolymer chain perpendicular and parallel to the interface in terms of chain length N and 
block size M, as well as the selectivity parameter x ■ The scaling relations differ for the case 
of weak and strong localization. In the strong localization limit a scaling relation for the lateral 
diffusion coefficient D\\ is also derived. We implement a dynamic off-lattice Monte - Carlo model 
to verify these scaling predictions. For chain lengths in a wide range (32 < N < 512) we find good 
agreement with the scaling predictions. 

PACS numbers: 82.35.Jk,83.80.Uv,07.05.Tp 
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I. INTRODUCTION 



The behavior of hydrophobic - polar copolymers (or HP - copolymer) at a selective 
interface (the interface which divides two immiscible liquids, say, water and oil, each of 
them being good for one type of monomers and bad for the other) is of great importance in 
the chemical physics of polymers. HP - copolymers are readily localized at such an interface 
because under a sufficiently large degree of selectivity the hydrophobic and polar parts of a 
copolymer chain try to stay on different sides of the interface due to the interplay between 
the entropy loss in the vicinity of the interface and the energy gain in the proper solvent. 
This reduces the interfacial tension between the immiscible liquids and motivates the main 
practical use of such copolymers as compatibilizersi. Early in the study of compatibilizers, 
the attention has been focused mainly on diblock copolymers due to their relatively simple 
structural characteristics2.i2.. These studies of diblock compatibilizers have been extended 
later by making use of the numerical self - consistent field approach and of Monte Carlo 
(MC) simulation o'^1^1^1 ''i^ . The seminal paper by Garel et. al.^ has attracted attention to the 
interface localization of random HP - copolymer and has brought about a great number of 
publications on this subjec t '^"I'^^i'^^i'^^i'^'^i'^^i'^^i^'^i^^i^^ . As a rule, the macromolecules considered 
in these studies pertain to random or strictly alternating copolymers whereas the impact of 
block size on the localization transition and on the properties of the interface have found 
comparatively little attention. In fact, to the best of our knowledge the only MC - simulation 
of multiblock HP - copolymers at a selective interface has been carried out by Balasz et al.— 
within a two-dimensional model where no treatment in terms of scaling has been suggested. 

In this paper we focus on the interesting special case of multiblock copolymers made 
of a regular sequence of HP - diblocks where the fractions of H - and P - segments are 
equal. The main objective is to bring together simple scaling arguments and the findings 
of MC - simulations. The principal idea of the scaling consideration in this case implies 
(see ref .^^1^^1^? ) that the problem at hand can be mapped onto the homopolymer adsorption 
problem whereby each HP - diblock (e.g. each part of consecutive blocks made of H- and 
P-monomers, respectively) acts as an effective segment. At the interface these HP - diblocks 
behave as "soft dipoles" so that their "polarization" minimizes the free energy. This phe- 
nomenon governs the HP - diblock effective attraction energy, and enables the mapping onto 
the problem of a homopolymer adsorption at a penetrable surface2i>2^. 
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In Sec. II we provide refined scaling considerations which are partially based on previous 
result si^i^iiS^. The results are presented in such a form that the chain and block length 
dependences are clearly separated and can be easily compared with the MC - findings. In 
Sec. Ill we report on the results of a comprehensive MC - simulation study and compare 
these results with the scaling predictions. We conclude with a brief summary and an outlook 
for further work. 

II. SCALING 
A. Static 

As mentioned before, the scaling treatment is based on the idea that a regular multiblock 
copolymer can be considered as a homopolymer where a single HP - diblock plays the 
role of an effective segment. All such diblocks try to keep their H - and P - segments in 
the corresponding preferable environment. This leads to the diblock "polarization" at the 
interface and to a free energy gain which produces an effective attraction energy E. We 
restrict ourselves here to the fully symmetric case. This means that the energy gain x 
(which is the selectivity parameter) for a P - monomer in its own (polar) environment is 
equal to the corresponding energy gain for a H - monomer, provided the latter stays in the 
hydrophobic environment. Also, the P - and H - parts of a diblock have an equal number 
of monomers. An estimate for the effective attraction energy per diblock of length 2M (M 
being the number of monomers of each species in the diblock) then reada^iiS^. 

E ~ -x'M\ (1) 

where we measure the energy in units of ksT, (i.e. x/^bT ^ x is dimensionless) where 
denotes the Boltzmann constant. 

Now consider a multiblock copolymer of length which consists of AT = N/2M diblocks. 
At the adsorption threshold, Xc, the attraction energy per chain scales as^ Vattr ~ \E\Af'^, 
where denotes the so called crossover exponent. It has been shown rigorously2£ that the 
relation (p = 1 — u (here u = 0.6 is the Flory exponent in 3D space) is correct for the case 
of a penetrable surface. At Xc the attractive energy, \4ttr = Kittr; should be of the order of 
the thermal energy, i.e. 

lE^lAf^-" ^ 1. (2) 
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Taking into account eq. (^J, the expression for the critical selectivity then yields 

Xe-M-^i+'^V^^-^i-^^/^. (3) 

Eq.Q refers to the chain as a whole at the localization threshold Xc- As the selectivity 
parameter x increases, each part of the chain (named "blob") containing g monomers accu- 
mulates an energy of the order of ksT. Note that the chain "neutrality" requires that g/M 
is always an even number. The blob length g may be determined from \E\{g / MY^'^ ^ 1, so 
that 

The localized chain as a whole can be considered as a string of such blobs. 

It is customary^'' to take the number of blobs, N/ g ~ ^^^/(i-i/) ^ natural 

scaling variable in the various scaling relations. For the sake of convenience, we take this 
variable in the following form 

Then the chain size perpendicular to the interface direction scales as 

Rg± = lN^G{yi) , (6) 

where I is the Kuhn segment and G{fi) is some scaling function. 

One should discriminate between the cases of weak and strong localization^^. In the weak 
localization regime the scaling argument r] is of the order of one (cf. eq.Q), and the scaling 
function Q{ri) has a power law behavior, i.e. 

, . I 1 , at 77 < 1 
Q{v) = { , (7) 

I ?7" , at 77 > 1 

where a = —2v/{l — v). This value of the exponent a has been found from the condition 
that at 77 > 1 the value of -R^j, is determined by the blob size alone, i.e. it does not depend 
on N. 

Similar reasoning can be used for the chain size in direction parallel to the interface 
namely 

= /iV^H(r/) , (8) 



where 7i{ri) is some scaling function. Again, the form of 'H{ri) can be determined by the 
condition that a.t rj > 1 the locahzed multiblock chain behaves as a two - dimensional self - 
avoiding string of blobs, i.e. -R[| ~ N"^, where z/2 = 0.75 is the Flory exponent in 2d space. 
As a result one finds 

{1 , at r? < 1 
' , (9) 

1]''^ , at 7] > 1 

where f3 = 2(z/2 — z/)/(l — z/). 

As the selectivity parameter x grows further, a characteristic point x = Xoo can be 
reached where Rg± and Rg\\ approach a plateau and do not change any further. This is the 
strong localization limit where the number of monomers in the blob becomes of the order 
of the block length, i.e. g ^ M. By making use of this condition in eq. we can write for 
the characteristic selectivity Xoo the following simple relation 

Xoo~M-\ (10) 

which corresponds to the well known Flory- Huggins limit of phase separation. Note that Xca 
depends only on the block size M and does not depend on the chain length A^. Given this 
value one can go back to eqs. (jHl - EI) and check the scaling in the strong localization limit. 
The resulting relations are 

Rg± = (11) 

= /MM — ~ M-^'^^-'^^AT^^ (12) 



which reflect the pancake geometry of a polymer localized at the surface. The relations (jlip 
and ()12|) can be expected because in this regime all P - and H - segments are predominantly in 
their preferred solvents. In addition, the scaling estimate of the polymer size perpendicular 
to the interface shows that the perpendicular extension of the chain is entirely determined 
by the size of the blob. 

The free energy gain in the localized state is proportional to the number of blobs22, i.e. 

where eq. (jH) has been used. In terms of r] eq. (|T!I|) reads 

Floe = r]'/^'-"^ at (14) 
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where r/oo = XooM(^+^)/2jY(i-i.)/2 ^ (Ar/M)^^-'')/^. At x > Xoo the free energy gain follows 
the strong localization law 

F,oc--xN , (15) 
which is mainly triggered by the energy gain. 

B. Dynamics 

After the localization on the selective interface, the multiblock copolymer can only diffuse 
in the two - dimensional space along the surface. Here we give the scaling estimate for the 
characteristic time r which is necessary for a chain to displace along the surface a distance 
equal to its own gyration radius. We define this time in terms of the chain length and 
the block length M in the strong localization limit. 

The diffusion coefficient of the whole chain along the surface can be estimated as 

/1\ 1 -Dblock 

where -Dbiock is the diffusion coefficient of one block, ( = (oM the corresponding friction 
coefficient and the number of diblocks Af = N/2M. The characteristic time of the 2 - 
dimensional chain displacement obeys 



T 



where the characteristic size of the block (which spreads in the 3 - dimensional space) 
-Rbiock ~ iM'^ {u = 0.6) whereas the 2 - dimensional Flory exponent is z/2 = 0.75. 

Now one can express all quantities in eq. ()17|) in terms of A^ and M. As a result one 
obtains 

^ Co ^'M-°-3 A^2.5_ (^g) 
It is clear that this result should hold for sufficiently long chains and blocks. 
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III. MC - SIMULATION 



A. Description of the model 

The off-lattice bead-spring model has been used previously for simulations of polymers 
both in the bulk2^>2S and near confining surface o^^i^-^i^^i^^i^^i^^ ; thus, we describe here the 
salient features only. Each polymer chain contains effective monomers connected by 
anharmonic springs described by the finitely extendible nonlinear elastic (FENE) potential. 

f/™E = -|i?^ln[l-^^^]. (19) 

Here ^ is the length of an effective bond, which can vary in between Crnm < i < C-max-, with 
^min = 0.4, Imax = 1 being the unit of length, and has the equilibrium value /q = 0.7, while 
R = imax — io = io — imin = 0.3, and the spring constant K is taken as K/ksT = 40. 
The nonbonded interactions between the effective monomers are described by the Morse 
potential 

Um = eM{exp[-2a(r - r„i„)] - 2exp[-a(r - r^m)]} , (20) 

where r is the distance between the beads, and the parameters are chosen as r^m = 0.8, 
ejii = 1, and a = 24. Owing to the large value of the latter constant, Uuij) decays 
to zero very rapidly for r > Tmini and is completely negligible for distances larger than 
unity. This choice of parameters is useful from a computational point of view, since it 
allows the use of a very efficient link-cell algorithn>2£. From a physical point of view, these 
potentials eqs. (fT!I|) . ()2()|1 make sense when one interprets the effective bonds as kind of Kuhn 
segments, comprising a number of chemical monomers along the chain, and thus the length 
unit imax = 1 corresponds physically rather to 1 nm than to the length of a covalent C — C 
bond (which would only be about I.SA). Since in the present study we are concerned with 
the localization of a copolymer at good solvent conditions, in eq. (j20p we retain the repulsive 
branch of the Morse potential only by setting Uuir) =0 for r > Tmin and shifting Uuir) 
up by Cm- 
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The interface potential is taken simply as a step function with amplitude Xi 

Uint{n,z) = { 

where the interface plane is fixed at z = 0, and ^"(n) 
distinguishes between P- and H- monomers. The energy gain of each chain segment is thus 
— X, provided it stays in its preferred solvent. 

Two typical snapshots for a chain with = 128 and M = 8 in the weak {x = 0.25) and 
strong [x = 10) localization limits are shown in FigQto illustrate the main features of the 
model. 

In each Monte Carlo update, a monomer is chosen at random and one attempts to 
displace it randomly by displacements Ax, Ay, Az chosen uniformly from the intervals 
—0.5 < Ax, Ay, Az < 0.5. The transition probability for such an attempted move is simply 
calculated from the total change AE of the potential energies defined in eqs. (fTI?l - l?H) as = 
exp{—AE/kBT). According to the standard Metropolis algorithm, the attempted move is 
accepted only if W exceeds a random number uniformly distributed between zero and unity. 
Since our potentials are constructed such that chains cannot intersect themselves in the 
course of random displacement of beads, one does need to check separately for entanglement 
restrictions. In the course of the simulation the starting configuration of the copolymer is 
relaxed for a period of time tq before measurements of the chain properties are carried out. 
We sample various static and dynamic quantities as the components of the gyration radius 
perpendicular, Rg±, and parallel, Rg\^, to the interface, the density distribution of the two 
kinds of monomers around the interface, internal energy, specific heat, diffusion coefficients, 
etc. We use periodic boundary conditions in the plane of the interface while there are rigid 
walls in the z - direction, where the simulation box extends from z = —16 to z = 16. The 
algorithm is reasonably fast, i.e. one performs ~ 0.5 x 10^ updates per CPU second on a 
2.8 GHz PC. Typically we studied chains with length 32 < < 512 and all measurements 
have been averaged over 2^^ samples. 



z>0 

(21) 

z <Q 

= ±1 denotes a "spin" variable which 
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B. Results 



The localization at the selective interface can be considered as a phase transition at least 
in the limits: ^ oo, M oo and x ^ with r]^ = ^^j\j{i-u)/2 j^{i+u)/2 ^ 
finite N, M and x the transition looks like a smooth crossover and one can define an order 
parameter (OP) in terms of the fractions of P - and H - monomers in the polar (at z > 0) 
and hydrophobic (at z < 0) semispace, respectively, i.e. 



Fig|21 shows the variation of this OP with the selectivity parameter x for ^ chain length 
= 128 and various lengths of the blocks M. Evidently, for sufficiently large M the OP 
reaches saturation at much lower values of x than for short blocks. It can also be seen that the 
chains with small M even in the localized state display a much lower degree of ordering than 
those with longer blocks since in the localized phase the number of "frustrated" monomers 
which stay in the wrong solvent increases with vanishing proximity to the interface. 

The transition location, Xc defined as the infiection point of the curves in Fig. Elcan be 
easily calculated. A more comprehensive information about the localization is revealed from 
the monomer density distribution histogram which is given in Fig. IHl While for x = 0.25 the 
copolymer is still not localized, even at x = 0.75 > Xc there is still a considerable amount of 
segments which stay in the wrong solvent and one needs a rather strong segregation x = 5.0 
in order to keep the segments entirely in their preferred environment. 

The localization at the interface is also marked by a maximum in the specific heat of the 
system - C„(x) = {NkhT^)~^ ((-^^) ~ (-^)^) - which we obtain from the fiuctuations of the 
internal energy of the copolymer E and show in Fig. |31 One can readily see from Fig. |3] 
that the specific heat curves collapse (despite significantly larger scatter in the measured 
data) onto different master curves for the two block sizes in concern. However, the heights of 
the C„-maxima at the localization transition appear rather insensitive with respect to chain 
length A^. Indeed, it has been shown earlier— that the specific heat critical exponent is zero, 
and there exist only weak logarithmic corrections to finite size scaling (C„ ~ [InNY^^^ for 
a Gaussian chain) . One can also infer from Fig. 0] a rapid decrease in the height of the 
Ct,-maximum with growing block size M which is to be expected bearing in mind that the 
chance for an energy fiuctuation rapidly diminish as less and less monomers cross into the 
"wrong" solvent with increasing M. 




(22) 
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1. Weak localization 

We turn now to the comparison of our data with the scahng prediction given in Sec. II. 
Figure displays the variation of the chain size perpendicular to the interface Rg± with 
scaling parameter rj. One can easily distinguish between the cases of weak and strong 
localization. In the first case the curves collapse nicely on a master curve in a range of r] 
which gets narrower with growing M. This proves the scaling law given by eq. (|7j) and 
is also in line with the trend given by eq. (fTUj) for the characteristic selectivity Xoo at the 
onset of the strong localization regime. The scaling function, eq. (|7j), is characterized by 
the exponent a = —3, which is indicated by the dashed line in Fig. Eland appears in good 
agreement with the simulation data. An interesting case which we observe in all our studies 
is shown in the inset of Fig. El- in the case of a diblock (i.e. when M = N/2) in the strong 
localization regime one finds Rg±{i]) / Rg±{0) > 1. Clearly, as the two blocks tend to stay 
away from the interface, thus maximizing entropy, they induce an elongation of the coil 
perpendicular to the interface in contrast to the behavior of multiblock copolymers. 

The set of graphs shown in Fig. El has been used to evaluate the critical Xc taken as 
the positions of the respective inflection points. The result of this evaluation is shown in 
Fig. El The insert shows the M - dependence of Xc at different A^. Evidently, the rescaling 
according to eq. Q leads to the master curve. The dashed line demonstrates the scaling 
prediction Xc ~ M~°'^ and is again in a good agreement with theory. Not surprisingly, at 
small M one observes deviations from the expected power-law behavior which demonstrate 
that scaling validity is achieved in the asymptotic limit only. 

Similar scaling representation has been done also for the parallel component of the chain 
size Rg\\ (see Fig. |7|). In the weak localization regime all data collapse again into a master 
curve. The scaling exponent (3 = 3/4 (noted by the dashed line, see eq.© is fairly close to 
the simulational results. 

2. Strong localization 

The area where Rg± and -R^n abandon the master curves in Figs. El HI and approach 
different plateau values at high t] corresponds to the strong localization regime. In this case 
the scaling should follow eqs. (fTTll and (fT^ . We have evaluated Rg±{oo) from the plateau 
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height values in Fig. El (dropping out the factor 1/Rg±{0)). The corresponding results are 
shown in Fig. |H1 The M - dependence of Rg±{oo)/l is given for different chain lengths N. 
The finite size analysis indicates that an extrapolation to — leads to the correct Flory 
scaling Rg^{oo) ~ M°'^ (see the dashed line in Fig. |H)). Fig. El makes it also apparent that 
the parallel component Rg\\{oo) follows to a good accuracy the law -Rg||(oo)// ~ M^^'^^N^-'''^ 
which comes from eq. (fT^ . Thus Fig. El demonstrates as well the scaling relation for the 
relaxation time r given by eq. ()17|) . The fit of our MC - results with the theoretical scaling 
prediction r ~ M'^'^N"^'^ is very good especially taking into account the relatively small 
values of the block length. 

IV. CONCLUSION AND OUTLOOK 

In the preset work we have studied the adsorption of regular HP copolymer at a pene- 
trable interface between two immiscible fluids. Specifically, we have focused our attention 
on the influence of the block size M on the adsoprtion transition as it is known that ef- 
ficient compatibilizing agents (i) reduce the surface tension and (ii) enhance the adhesion 
between the surfaces that result upon cooling. While reduced surface tension is achieved by 
possibly higher surface coverage by the polymer, for enhanced adhesion between the com- 
ponents each block should extend and penetrate a significant distance into the respective 
compatible melt. Qualitatively it is clear that shorter blocks will provide a more efficient 
reduction of the surface tension whereas longer blocks would favor better adhesion between 
the segregated systems. In order to get more quantitative insight into the effectiveness of 
multiblock copolymers we develop a simple scaling treatment which is then compared against 
simulational results. 

Within the framework of our simple scaling approach we arrive at several important con- 
clusions characterizing the impact of block length M on the behavior of regular multiblock 
copolymers at a fiuid-fiuid interface: 

• The critical selectivity decreases with growing block length as Xc ~ M~'^^+'^^/^ while 
the crossover selectivity to the strong localization regime vanishes as Xoo ~ M~^. 

• The size of the copolymer varies in the weak localization regime as Rg± ~ 
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• In the case of strong segregation we obtain respectively Rg± ~ M" and Rg\\ ~ 

• The typical relaxation time in the case of strong localization varies as r ~ M'^^'^~'^^^ 
with block length. 

Our computer experiments appear to confirm nicely these predictions. 

Concluding, one should point out that this study has been focused on one aspect of the 
localization of copolymers at penetrable interfaces. A number of related topics such as the 
case of many chains in the system, the presence of attractive interactions between similar 
segments, or the kinetics of localization could provide further valuable insights on the prob- 
lem. Another interesting aspect comes when random - or specifically sequenced copolymers 
are considered rather than regular multiblock copolymers. In these cases additional infor- 
mation is stored in the sequence of the H and P monomers. The relevant information can be 
read by the interfacial behavior of such polymers. Indeed, the localization will then depend 
on the distribution of the H and P monomers along the chain. For special sequences we 
expect to relate their behavior to some typical proteins^. We hope to be able to report on 
some of these topics in the near future. 
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FIG. 2: The order parameter for chain length = 128 vs selectivity x for different values of the 
block length M. 
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FIG. 4: Specific heat (per monomer)of a copolymer chain with length 32 < < 512 and two block 
sizes M = 2, 4 at the localization transition. 
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FIG. 5: Variation of the chain size perpendicular to the interface Rgj_ with the scahng variable 
rj for block lengths M = 1,2,4,8 and 16 (inset). In all cases the dashed line indicates the scaling 
prediction. 
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FIG. 6: Scaling plot of the critical selectivity Xc in terms of block length M and chain length A'^. 
The dashed line indicates the scaling prediction. The inset shows Xc vs M without scaling in N. 
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FIG. 7: Variation of i2g|| with the scahng variable rj for chain lengths 32 < < 512 and several 
block lengths, as indicated. The onset of the strong localization regime is manifested by the decay 
of the master curves into different curves for each N. In all cases the dashed line indicates the 
scaling prediction. 
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FIG. 8: The perpendicular component Rg±{oo)/l as a function of M at the different A^. The 
dashed hne corresponds to the scahng prediction. The insert shows the extrapolation for — > 
as leading to the theoretical value of the slope v = 0.6. 
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FIG. 9: The upper curve represents the rescaling of the characteristic time r according to eq. (|17() . 
The lower curve shows that the scaling for Rg^^{oo) follows ea.(|12[). The dashed lines indicate the 
scaling predictions. 
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